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We perform a Monte-Carlo study of the lattice two-dimensional gauged XY-model. Our results 
confirm the strong-coupling expansion arguments that for sufficiently small values of the spin-spin 
coupling the "gauge symmetry breaking" terms decouple and the long-distance physics is that of the 
unbroken pure gauge theory. We find no evidence for the existence, conjectured earlier, of massless 
states near a critical value of the spin-spin coupling. We comment on recent remarks in the literature 
on the use of gauged XY-models in proposed constructions of chiral lattice gauge theories. 



I. MOTIVATION AND SUMMARY 

The gauged two-dimensional lattice XY-model has 
been studied for quite some time. Its simplicity allows for 
analytic studies via lattice dualities and strong-coupling 
expansions [1 as well as for numerical analysis, in either 
a Hamiltonian [5] or Euclidean (see, e.g., [3]) formulation. 
The model exhibits many phenomena found in more com- 
plicated higher-dimensional theories, notably the pres- 
ence of both confining and Higgs phases [TJ HJ HJ [5] . 

The motivation for this short study stems from our in- 
terest in using the gauged XY-model, as well as its non- 
abelian higher-dimensional analogues [4j |6j [7] , in lattice 
constructions of chiral gauge theories. The lattice for- 
mulation of chiral gauge theories is an outstanding prob- 
lem with no practical solution yet, despite much recent 
progress; for reviews, see [8]. While not our topic here, we 
note that the recent constructions of [HI [TU] have some at- 
tractive features (for an earlier proposal of similar flavor, 
see [H]). Most importantly, they may offer a way around 
the difficult and unsolved problem of the explicit nonper- 
turbative construction of the fermion measure for general 
chiral lattice gauge theories. Ref. [9 aims to achieve this 
by combining older ideas [T2] to use non-gauge strong 
dynamics to decouple the mirror fermions in a vectorlike 
theory with the recently discovered exact lattice chirality 
of the Ncubcrgcr-Dirac operator. 

A concrete two-dimensional realization of the proposal 
[3] has the gauged XY-model as an essential part and 
was studied in various limits in [10 . Recently, |13j con- 
sidered the perturbative spectrum of the models of [5] at 
nonzero gauge coupling. It was claimed there that the 
gauge boson is always massive and, therefore, the con- 
struction of [U [10] was argued to be irrelevant to the 
study of unbroken chiral gauge theories, on account of 
this fact alone. 

This brings us to the main topic of this paper. The 
arguments of [13] do not invoke either the fermions or 
the strong mirror dynamics (admittedly, complicated and 
not yet completely understood), but concern only the 
spectrum of the gauged XY-model, arguing that it only 
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has a phase with a massive gauge boson. This conclusion 
contradicts strong-coupling expansion arguments for the 
decoupling of the "gauge-breaking" terms for small values 
of the spin-spin coupling. As these arguments have been 
made many times in the past [H [5] , we will not repeat 
them at length here. 

Instead, given that confusion around the issue appears 
to still persist, we use Monte-Carlo simulations to argue 
the same point. These methods allow us to also study 
regions beyond the reach of strong-coupling expansions 
or perturbation theory and to observe in detail how the 
infrared physics changes smoothly as the XY-model spin- 
spin coupling k is varied. We show that for k < k c , with 
k c = 0(1), the long-distance physics of the gauged XY- 
model is that of the unbroken pure gauge theory. We 
numerically demonstrate that, in accord with the strong- 
coupling expansion, for subcritical k the leading effect of 
the spin-spin coupling can be incorporated into a shift of 
the gauge coupling. 

We also study the spectrum in coupling regimes in- 
accessible by controlled analytic methods. The lattice 
Hamiltonian of the gauged XY-model was studied long 
ago via a strong-coupling expansion in the gauge cou- 
pling, valid for arbitrary values of k, with a subsequent 
Pade resummation to small gauge coupling [3J. A dip 
in the spectrum of charged particle pairs near a critical 
value of k of order unity was found, and it was conjec- 
tured that it may indicate the presence of massless states. 
We look for such states by measuring the susceptibilities 
of the relevant operators. We find a slight increase of 
the susceptibility around similar values of k, consistent 
with the results of the Pade resummation of the strong- 
coupling scries. However, we find no finite-size scaling 
of the susceptibility and hence no evidence for massless 
states in the intermediate "critical" regime. 

Finally, as far as lattice formulations of chiral gauge 
theories are concerned, our results address the objection 
of [T3] to the proposal [S]. Much remains to be done to 
see whether constructions along the lines of 11J, |9J fulfill 
their designed goal — to give rise to long-distance gauge 
theories with chiral fermions in complex representations. 
While we hope to report more on at least some of the 
outstanding issues in the future, our main point here is 
that the reason for their failure would have to be more 
subtle. 



II. THE MODEL 



The Euclidean action of the gauged XY-model is: 



S 




(j)* x U Xifl (f) x+(l +h.c, (1) 



where U x ^ = e lAx ^ is the link group element, Y[ U is the 

plaq 

usual plaquctte Wilson action, and (f> x — e lr,m is a unitary 
Higgs field; x denotes lattice sites on a two-dimensional 
square lattice and fx is a unit lattice vector in the fj, = 1, 2 
direction. Here A x ^ and rj x are angular variables defined 
on links and sites, respectively; throughout the paper the 
lattice spacing is set to unity. 

The partition function of the model is defined via a 
lattice path integral over all angular variables. The ac- 
tion can be more succinctly expressed (and the simula- 
tion simplified) if one notes that upon a change of the 
link variables the dependence of the action on r\ x can be 
completely eliminated, leaving us with a simpler theory 
to study: 



S = Y,\ I3cosF x 



n cos A 
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where F x = A x+lt2 - A x<2 - A x+2: i + A x>1 is the lattice 
gauge field strength. In the naive continuum limit, taking 



the action (|2|) becomes 
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and describes a massive gauge boson of mass = 
g 2 K — a conclusion which holds for sufficiently large k. 
For small k, however, the perturbative expansion around 
the naive continuum limit is not a good guide to the spec- 
trum of the theory, because the naive continuum limit 
misses two important points: the angular nature of the 
variables and the fact that at small n the quantum fluc- 
tuations of the longitudinal component of the gauge field 
A x ^ are not suppressed. For n <C 1, the 0(A i ) terms 
lead to strong interactions between the longitudinal com- 
ponent of A^ — the fluctuations of these modes are not 
suppressed by the gauge invariant kinetic term and for 
small k the spin-spin (or "mass" ) term does not suppress 
them either. In contrast, for k> 1, these modes' fluctu- 
ations are suppressed by the "mass" term in ^ and the 
{A^) A interactions are irrelevant. 

For small k, therefore, a strong coupling expansion in 
k is more appropriate than the perturbative expansion — 
the longitudinal components of the gauge field strongly 
fluctuate on the scale of the lattice spacing, are thus 
heavy, and can be integrated out. The small- k expansion 
is most conveniently performed in the original gauge in- 
variant formulation in terms of the variables A x ^ and r\ x 



of Q. To leading order in k, the effect on the infrared 
physics of the gauge field of the rapid fluctuations of the 
Higgs field <j) x is to renormalize the gauge coupling [HE]. 
The long-distance — compared to the correlation length 
of the Higgs field excitations, which is smaller than the 
lattice spacing for k < k c — effective action is of the form: 



*S'k< 
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0- 
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where dots denote higher-dimensional gauge invariant 
terms. Computing the coefficient of the leading-order 
shift of the bare gauge coupling in Q , as well as those of 
the higher-order terms (leading to larger plaquettes and 
omitted in Q), is trivially done by using the small-K ex- 
pansion of the partition function defined by ([TJ , and by 
integrating over the fluctuations of the compact scalar 
field (j) x . 

Equivalently, we show in the next Section that the co- 
efficient of k 4 in Q can be read off the results of our 
simulation (clearly, this is numerically possible only if (3 
is not too large). The Monte-Carlo methods also allow 
us to see how the infrared physics changes as a function 
of k and to study the spectrum of various operators in 
regimes where neither the strong-coupling expansion nor 
perturbation theory apply. 



III. THE MONTE-CARLO STUDY 

Our Monte-Carlo study focuses on the evaluation of 
several quantities: the Wilson loop, the Polyakov loop, 
and several zero-momentum correlators, or susceptibili- 
ties. We study their variation with k, /3, and the lattice 
size N, our primary interest being the k dependence for 
small gauge coupling. 

This choice of observables is motivated by the nature 
of two-dimensional physics: we remind the reader that 
the two-dimensional pure U(l) gauge theory (which we 
argue describes the infrared of the gauged XY-model for 
subcritical k) has no local propagating degrees of free- 
dom. Nonetheless, it exhibits nontrivial dynamics mani- 
festing itself in a nonzero expectation value of the Wilson 
loop, with an area law behavior for all values of the gauge 
coupling. Similarly, the Polyakov loop expectation value 
vanishes in the pure gauge theory for any coupling and 
temperature. On the other hand, the gauged XY-model 
with nonzero sufficiently large k has propagating local 
degrees of freedom, and information on their spectrum 
can be inferred from correlators of local operators. 

All our simulations are performed via the Metropo- 
lis algorithm, using the partition function in "unitary 
gauge" defined by the action The autocorrelation 
time for the various observables is monitored. The largest 
autocorrelation time that of the Polyakov loop, varying 
between about ten and a few hundred of lattice updates, 
depending on f3 and k — is used to estimate the errors. 
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A. The Wilson and Polyakov loops 



String tension, N=8 



We denote by Wij the expectation value of the Wilson 
loop operator for a rectangular loop of size i X j. In the 
pure gauge theory, k = 0, the Wilson loop expectation 
value can be exactly calculated. The string tension a of 
the pure U(l) gauge theory, which confines for all values 
of the gauge coupling, see e.g. [13], is easily shown to be, 
in the infinite volume limit: 



i x j 



= - In 



HP) 



(5) 



where I±, Iq are Bessel functions. At strong coupling, the 
string tension is a = hxAg + (D(l/g 2 ), while at weak 
coupling one instead obtains a ~ g 2 . 

Finite-volume corrections to the string tension will 
sometimes be significant in our simulations. In order to 
disentangle these effects, we have computed the Wilson 
loop expectation value in the pure gauge theory in finite 
volume, using the dual representation (see, e.g., the sec- 
ond reference in 1 ) of the partition function. We thus 
obtain for IF,-.,- in finite volume: 




FIG. 1: The dependence of |lnWjj|/(y) on k for i = j = 
1,2,3,4, for an 8 2 lattice with (3 = 4. The figure clearly 
shows a transition to an area law at k < 0.5, with string 
tension given by the expression ([5| for the pure gauge theory, 
cr(4) « .147. A perimeter law is approached at large k. 



N 2 



(6) 



where A — ij is the area of the loop, TV 2 is the lattice 
volume, and the sums are over all integers. It is easy to 
check a few simple limits: as N — > oo |6]) reproduces the 
result ([5) for the string tension; also, it correctly yields 
Wnn — 1- We use ^ to show that for subcritical n 
our measured Wilson loop expectation values reduce to 
the ones of the pure gauge theory also in cases where 
finite-size effects are numerically significant. 

Turning on re 7^ 0, for large k we expect from Eqn. ([3| 
to obtain a perimeter law for the Wilson loop. For 
small re, from Eqn. Q we should expect an area law 
instead. To avoid confusion, we stress that the theory 
does have charged matter excitations — the quanta of 
the Higgs field <p x — which can screen external charges 
and lead to perimeter law for the Wilson loop. In the 
small-re limit, however, the charged excitations are heav- 
ier than the UV cutoff. While for any charged fields' 
mass a sufficiently long "string" will eventually find it 
energetically favourable to break, we can imagine taking 
the continuum and large- volume limits so that the break- 
ing does not occur even for the longest possible string, 
i.e., aN/2 < nig, where N is the lattice size and rog is the 
lowest energy of a pair of charged particles, much higher 
than the UV cutoff (mg ~ re -1 in the strong-coupling 
expansion [2]). 

We begin the discussion of our Monte-Carlo re- 
sults by showing the dependence of the "string ten- 



sion" aij——lnWij/(ij) on re for Wilson loops with 
j=j'=l,2,3,4. On Fig. [Tlwe show a for an 8 2 lattice 
and (3 = 4, while on Fig. |2| we show a on a 16 2 lattice 
for the same j3 = 4. For these loop sizes and values of (3 



String tension, N=16 




FIG. 2: Same as Fig. [T] but for a 16 2 lattice with f3 = 4. 



the finite volume effects are small. For small re < 0(1), 
we see that in both cases, the pure gauge theory value 
of the string tension er(4) ~ .147, as given by Eqn. (|5j, 
is reproduced. For large values of re ^ 1, on the other 
hand, we see a gradual transition to a perimeter law for 
the Wilson loop, with the "string tensions" for re = 4, 5, 8 
scaling as inverse powers of \J~A~ (re = 8 comes closest to 
pure perimeter law). 

The strong-coupling expansion Eqn. Q, predicts that 
for sufficiently small re the Wilson loop expectation value 
should be given by the same equation as in the pure gauge 
theory, Eqn. , but with (3 — > (3 + re 4 /8, so long as the 
loop is larger than the correlation length of the massive 
charged excitations. The small effect due to the C(re 4 ) 
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Polyakov line, N=16 




FIG. 3: The fit of A/3(k) vs. k 4 /8, see Eqn. Q for the Wilson 
loop expectation value. Fit uses an 8 2 lattice, 1x1,2x2 and 
3x3 loops. 



FIG. 5: Expectation value of the Polyakov loop on a 16 2 
lattice as a function of k. A transition to a confining phase 
for «: < 1 is clearly visible. 



Wilson loop expectation value, N=8 



l 2 , 2 2 , 3 2 , and 4 2 , as function of k. The prediction for 
Wu of the pure gauge theory, including finite-size effects, 
Eqn. Q, is given in the caption of Fig. [4] for k < 0.7, 
the agreement with the numerical results is evident. 

Finally, we study the expectation value of the Polyakov 
loop operator around the /x-th compact direction, P M ; we 
recall that its vanishing signals the onset of a confining 
phase where the free energy of a isolated quark is infinite. 
On Fig. [5j we show (the loops in the two directions 
are identical, up to numerical errors) measured on a 16 2 
lattice, for a range of values of /3, as a function of n. As 
with the Wilson loop, the fact that the Polyakov loop 
achieves the value appropriate to the pure gauge theory 
(zero) at k < 1 is clearly visible on the figure. 



FIG. 4: Expectation values of Wilson loops of sizes 
1 2 ,2 2 ,3 2 ,4 2 , for ft = 20, for an 8 2 lattice as a function 
of k. For k < 0.5 the results are in agreement with the 
finite-size formula Q for the pure gauge theory, which gives: 
W11 « .975, W22 » .908, W 33 « .819, and W44 » .732. 



shift of the bare gauge coupling can not be seen for the 
relatively large values of (3 on Figs. [TJ [2| To test this 
prediction and see for what values of k it breaks down, 
we have measured Wu, W22, and W33 as a function of 
/3, near j3 ~ 1.6, for a number of values of k, ranging 
from —1 to 1 with spacing 5 k = 0.1. Then, we determine 
the function (3(k) such that Wu is constant. On Fig. [3] 
we plot A/3 = P(k) — (3 vs. the small-/t prediction of 
the strong-coupling expansion, A(3 str . coup i. ~ k 4 /8. The 
agreement for k < 0.5 is quite persuasive. 

We present one more plot of the Wilson loop expec- 
tation value, this time for the larger value of (3 = 20 
and an 8 2 lattice, on Fig. H where we show Wu of sizes 



B. Zero-momentum Green's functions 

In this Section, we study the following zero-momentum 
correlators (susceptibilities) : 

X 

X^u = ^(sin^^sinAo^) , (7) 

X 

where (■■■) c denotes a connected correlator; in this Sec- 
tion, we take the indices to be \x, v — 0, 1. 

We use the susceptibilities to look for massless states 
in the relevant scalar (S) and vector (V) channel via 
a large- volume scaling. In a Hamiltonian formulation, 
the operator Ke((f)*U x j0 x , j), equal to cosA^i in uni- 
tary gauge, creates pairs of oppositely charged particles, 
while lrn(tf>* x U x i4> x+ i) = sinA^i creates (massive) vec- 
tor states. The spectrum of the lattice gauged XY-model 
Hamiltonian was studied in ref. [5] via an expansion in 
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Susceptibilities beta=4 
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ond, the 00 and 11 components are identical due to the 
same rotational symmetry of the Euclidean action; 
hence, we only show the 00 component. Third, it is clear 
from the figure that there is no large-volume scaling of 
the vector susceptibility for any value of k, and hence no 
indication that massless states are present; we have also 
checked that this lack of scaling persists for larger values 
of /3. Finally, as explained in Section [TTJ for k — > oo per- 
turbation theory is valid and the susceptibility ^ can 
be approximated by: 



Xoo 



' lim 



9oo 



PaPa 
g 2 K 



g 2 K 



~, (8) 

K 



FIG. 6: The values of the vector susceptibilities Xoo (upper 
set of curves) and Xoi (vanishing lower curves) as functions 
of k for lattices of size 4 2 , 8 2 , 16 2 . The 1/k behavior expected 
in perturbation theory, see Eqn. valid for large k, is also 
indicated. 



Susceptibilities beta=4 



<cos-cos> N=4 > 
<cos-cos> N=8 t 

<cos-cos> N=1 6 : 
<cos-cos> N=4 i 
<cos-cos> N=8 t 

<cos-cos> N=16 L 




FIG. 7: The two eigenvalues of the scalar susceptibility matrix 
xL as functions of k for lattices of size 4 2 ,8 2 , 16 2 . Clearly, 
no large- N scaling is observed near the transition. 



1/g 2 , valid for all values of n. Subsequently, a Pade re- 
summation was used to continue the series to g 2 — > 0. 
The [4, 4] Pade approximant of the scalar correlator com- 
puted there showed a tendency for scalar (S) states to 
become light for k near n c (see Figs. 8, 9 of ref. [5]). 
While ref. [2J admitted that it is not clear that the pro- 
cedure used to arrive at this result is to be trusted, it 
conjectured that there may be massless states near k c . 

We first show our results for the vector, (sin sin A v ) , 
susceptibilities (|7| on Fig. [6j for lattices of various sizes 
and for /3 = 4. A few remarks are in order. First, the off- 
diagonal (fif = 01) components of the vector correlators 
vanish due to the Aq — > A\, A\ — > —Aq symmetry. See- 



where we used the perturbative expansion of the ac- 
tion, Eqn. ([3]), to calculate the zero-momentum two-point 
function of the massive gauge boson (the overall factor 
of g 2 is due to the field normalization and the metric 
is Euclidean). Comparing the results of the simulation 
with the 1/k curve, also drawn on Fig. [6] we see that the 
expected perturbative behavior ^ is well reproduced by 
our Monte-Carlo results for large k, and that the Green's 
function is drastically different for small values of k. 

The measurements of the scalar susceptibilities, 
(cos A^ cos A v ), are shown on Fig. [7] The states cre- 
ated by co$,A x ^ and cosA^i mix, as the Z4 symmetry 
does not forbid them to (in the Hamiltonian formulation, 
the mixing is between states created by the operators 
Re(0*J7 a , i</) x+ i) and cos tt x ) . We show the two eigenval- 
ues of the susceptibility matrix on Fig. [7] Clearly, while 
an increase of the susceptibility (and hence the corre- 
lation length) around the critical value of k is visible 
on Fig. [7] — consistent with the Pade resummation of the 
1/g 2 expansion [2] — the absence of finite-size scaling in- 
dicates that the two-particle scalar states remain heavy, 
with a correlation length of order the lattice spacing. 

While our main interest in this Section was the sus- 
ceptibility (the zero-momentum Green's function) and 
its large-volume scaling, we have also measured position- 
dependent Green's functions. The results reveal that, for 
k < 1, the fields are essentially uncorrelated at neigh- 
boring lattice sites, and are thus heavier than the UV 
cutoff; we note that this is also the result obtained in [2], 
We do not show these data as they are not particularly 
illuminating; moreover, even if one takes the value of the 
zero-momentum Green's functions for k < 0(1) as in- 
dicative of the square of the Compton wavelength of the 
state (this, in fact, is a gross overestimate, as the true en- 
ergies of the states go to infinity as n — > 0) , a quick look 
at Figs. [7J [6] shows that the maximum Compton wave- 
length is about a/O^ w .84 in units of the lattice spacing. 
This is consistent with the strong-coupling expansion in 
k, as described at the end of Section [TTJ, as well as with 
the results of ref. [2], which showed that both V and S 
states decouple at k < k c . 
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IV. CONCLUSIONS FROM THE 
MONTE-CARLO STUDY 

We performed a detailed study of a number of observ- 
ables in the gauged XY-model, for various values of (3 
and k. Here is a summary of the results of this paper: 

1. We showed that for k < k c the string tension pre- 
cisely agrees with the exactly known pure gauge 
theory expression, Eqn. ([5]), for a number of lat- 
tice sizes and values of f3, see Figs, [I] [2] [4] This 
even includes observable finite size effects in Fig. [4] 
which correctly reproduce the finite-volume pure 
gauge theory result of Eqn. Q. 

2. In Fig. [3j we observed even the small C(k 4 ) correc- 
tion to the bare gauge coupling of Eqn. Q from the 
leading-order strong-coupling (small- /c) expansion. 

3. We showed that the Polyakov loop expectation 
value vanishes for k in the same range, k < n c , 
where the string tension is that of the pure gauge 
theory, see Fig. [5] 

4. Our results for the scalar, Fig.[7J and vector, Fig. [6] 
susceptibilities ^ revealed no evidence for mass- 
less states in the spectrum in the intermediate 



regime of k ~ k c . 

5. Furthermore, the vector susceptibility was shown, 
see Fig. [6j to match that expected of a massive 
gauge boson for large values of k > 1, precisely 
where perturbation theory is expected to be valid, 
see Section [nj The same zero-momentum propa- 
gator is, however, drastically different for k < k c , 
indicating that the correspondent state is heavier 
than the inverse lattice spacing and decouples from 
the infrared physics. 

These results show that we have achieved our main ob- 
jective: to provide numerical evidence for the fact that 
at small n < k c — 0(1) the "gauge symmetry breaking" 
decouples and that the infrared physics of the theory is 
that of the pure gauge theory, with a small shift of the 
bare gauge coupling. 
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